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ABSTRACT 

We have obtained time scries optical spectra of the cataclysmic variable PQ And in quiescence. 
The spectra show a white dwarf continuum with narrow Balmer emission superimposed over 
strong Balmer absorption. The emission lines have blue and red components whose strength 
changes with time. An analysis of the Ha emission line implies a short orbital period below the 
period gap. Given its lack of accretion disk features, its large and infrequent outbursts, and an 
orbital period below the period gap, PQ And is probably a low accretion rate object similar to 
WZ Sge. In addition, white dwarf model fits imply that PQ And is an excellent ZZ Cet candidate. 

Subject headings: binaries: close — stars: dwarf novae — stars: individual (PQ And) 



1. Introduction 

PQ And was discovered by McAdam on 21 
March 1988 (Hurst et al. 1988a) at a visual magni- 
tude of 10. Examination of the Palomar Sky Sur- 
vey plates showed a precursor object with a blue 
magnitude between 18-19 and a red magnitude 
of about 20 (Hurst et al. 1988b; Hurst & Young 
1988). Within 19 days the light curve had declined 
by 2 magnitudes. The large outburst amplitude 
and rapid decline led to the initial classification 
of a classical nova but spectra taken 3.5 months 
after the outburst by Wade & Hamilton (1988) 
lacked the nebular features typically observed in 
classical novae. Instead, the spectra had strong 
Balmer emission surrounded by broad absorption 
similar to the quiesent spectra of the dwarf nova 
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WZ Sge. The strong absorption indicated that the 
energy distribution was primarily from the white 
dwarf (WD) and implied an extremely low accre- 
tion rate. Recently, it has been shown that sys- 
tems with this type of spectra are good candidates 
for containing non-radially pulsating DA WDs or 
ZZ Cet stars (Woudt & Warner 2004; Warner & 
Woudt 2004). 

WZ Sge is a dwarf nova whose outbursts occur 
on timescales of tens of years and with amplitudes 
of 7-8 magnitudes. Cataclysmic variables (CVs) 
of this type are also known as TOADS or Tremen- 
dous Outburst Amplitude Dwarf novae (Howell, 
Szkody, & Cannizzo 1995). These characteristics 
are thought to be due to the very close binary 
separation and the extremely low accretion rate 
in these systems. TOAD orbital periods are all 
below the period gap and are typically of order 90 
minutes. The mass transfer rates implied by the 
models are only 10 _11 M Q yr _1 (Howell, Szkody, 
& Cannizzo 1995). 

Given PQ And's similarities to a dwarf nova, 
Richter (1990) searched archival plates for evi- 
dence of previous outbursts. Two other outbursts 
with maximum magnitudes of ~ 11 were found on 
23 August 1938 and 7 March 1967. With decades 
long superoutburst timescales, an outburst ampli- 
tude of ~ 9 mag, and a WZ Sge like quiescent 



spectrum, PQ And clearly has TOAD character- 
istics. In this paper we show that the orbital pe- 
riod appears to be below the period gap and that 
the system likely contains a WD in the instability 
strip. 

2. Observations 

We obtained 10 spectra of PQ And on the 6.5m 
Multiple Mirror Telescope (MMT) on 16 Septem- 
ber 2003. We used the blue channel CCD spectro- 
graph with a 500 grooves per millimeter grating 
and a 1" wide slit for a spectral resolution of 3. 6 A. 
Spectral coverage extends from 4100A to 7200A. 
The exposure time was set to 480 seconds in order 
to resolve an orbital period as low as 80 minutes. 
The seeing ranged from 1'.'0-1'.'5 during the night. 
Biases and flat fields were obtained at the begin- 
ning and end of each night and a HeNeAr compar- 
ison lamp spectrum was obtained after each PQ 
And observation for optimum wavelength calibra- 
tion. Standard stars were also observed through- 
out the night for flux calibration. The spectra were 
reduced using standard IRAF 1 routines. 

The combined spectrum is shown in Figure l 2 . 
The combined spectrum has narrow emission lines 
flanked by broad absorption similar to that de- 
scribed by Wade & Hamilton (1988). The equiv- 
alent widths, fluxes, and FWZI values of the ob- 
served emission lines in the combined spectrum 
are given in Table 1. All the emission lines have 
a red and a blue component whose strength and 
peak vary with time. The evolution of these com- 
ponents can best be seen in Ha and is shown in 
Figure 2. 

3. Emission line analysis 

Radial velocities from the line wings were mea- 
sured by convolving the line profiles with a pair of 
Gaussian bandpasses (see Shatter 1985; Szkody et 
al. 2002, for details) whose width is set to the res- 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 

2 Since the extinction toward PQ And provided by the 
Schlegel, Finkbeiner, & Davis (1998) maps is extremely 
low, E(S — V) = 0.06, none of the spectra shown in this 
paper have been dereddened nor have they been Dopplcr 
corrected. 



olution and whose separation is varied. For each 
separation, the velocities are fit to a sine function 
of the form 

v = 7 - ifsin[27r(</> - <j) )] (1) 

where 7 is the systemic velocity, K is the semi- 
amplitude, and <p is the offset from superior con- 
junction of the WD. The best solution is chosen 
as the one producing the smallest errors in the in- 
dividual parameters and the total a of the fit. 

We fit both the Ha and H/3 lines with this 
method but unfortunately not all the observed line 
profiles were useful. Two Ha and one H/3 profiles 
were discarded due to cosmic rays or noise at a par- 
ticular separation. Even with one less data point 
the Ha fit was significantly better than that deter- 
mined for H/3 thus we adopt the Ha separation and 
period as the best solution and force H/3 to match. 
The best fit period is 1.7±0.1 hours where this is 
the formal uncertainty from the analysis. The true 
uncertainty in the period is greater since the de- 
rived period is longer than the sampling time. Fig- 
ure 3 shows the radial velocity data for Ha and H/3 
fits. The radial velocity data is provided in Table 
2. The orbital parameters of the lowest total a fit 
are given in Table 3 along with other parameters 
of PQ And. 

To quantify the component changes, we calcu- 
lated the ratio of the blue (violet) to red compo- 
nents (V/R) for each spectrum that shows Ha. 
The V/R ratios were calculated with the IRAF 
program splot with two Gaussians which had the 
same FWHM and were centered at 6554Aand 
6574A. In each spectrum the Gaussians were only 
allowed a single wavelength shift provided that the 
separation was always 20A. The results are shown 
in Figure 4. 

Superior conjunction of the WD is defined as 
the red to blue crossing point of the emission lines. 
In Figure 4 we have fit the V/R data with a sine 
wave of the form: 

V/R = A + Bsin[27r(t - t /P O rb)] (2) 

to determine the crossing point. With P or b = 1.7 
hours the best fit gives a superior conjunction date 
of HJD 2,452,898.902. The other model parame- 
ters are provided in Table 3. 
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3.1. Differential photometry 

As an added check on the derived period we 
utilized differential photometry of PQ And from a 
CV orbital period survey taken in early 1999. The 
data were obtained at the Steward Observatory 
1.6m Kuiper Telescope on Mount Bigalow with the 
2Kx2K back-illuminated CCD Imager. The CCD 
was set at 2x2 binning giving an image scale of 
0.3" /pixel. The seeing varied during the run be- 
tween 1" and 2". The observations were carried 
out in the / band where the accretion disk and 
secondary star are typically of comparable bright- 
ness. The exposure time in the PQ And field was 
600 seconds. Biases and flats were taken at the be- 
ginning and end of each night and the data were 
reduced in IRAF. The only night that we spent 
sufficient time to confirm a 1.7 hour period was 
19 January 1999. On that night PQ And was ob- 
served 18 times over 3 hours. A differential light 
curve of the results is shown in Figure 5. The ob- 
served changes are suggestive of a variation similar 
to our derived spectroscopic period. The normal- 
ized pcriodogram (Scargle 1982; Home & Baliunas 
1986) of the photometry gives a longer period of 
2.1 hours but at a low confidence level of only 50%. 

We emphasize that these arc not definitive or- 
bital period solutions since the derived periods are 
similar in length to the total time that PQ And 
was observed. More data are necessary to reach 
a more definitive orbital period for this object. 
However, the fact that the radial velocities and 
the changes in the blue and red components show 
sine wave like behavior on the timescale of our 
observation runs imply that the period must be 
short. A firm upper limit on the orbital period 
can be established since the accretion disk would 
overwhelm the WD at periods > 3 hours and there 
are no known dwarf novae in the gap. Thus for the 
remainder of the paper we will assume an orbital 
period range between 1.7 and 2.1 hours. 

4. White dwarf model fitting 

The WD effective temperature and gravity were 
determined by comparing the observed spectrum 
to a grid of synthetic spectra computed using the 
PHOENIX model atmosphere code. Details of the 
modeling can be found in Barman et al. (2000). 
The grid includes models with T eff = 9,000K - 
15,000K (in steps of 1000K), log(.g) = 6.5, 7.0, 7.5, 



8.0, and 10 _1 , 10~ 2 and solar metallicities. After 
masking the Balmer emission lines due to the ac- 
cretion disk, T e f f was first determined using least- 
squares minimization across the full observed spec- 
tral range. Fixing T e ff, the surface gravity was 
found by fitting only the wings of H(5, Hy, and H/3 
since these were the only features available that 
arc reasonably sensitive to gravity The models 
with sub-solar metallicities gave only marginally 
better fits; however, models with solar metallic- 
ities indicate that several metal absorption lines 
should have been detected but clearly were not 
observed. Consequently, the metallicity is likely 
less than 10 _1 solar. Overall, varying metallicity 
had little effect on the final T e ff and log(g) values. 
The final best-fitting model had T eff = 12,000K 
± 1000K and log(g) = 7.7 ± 0.3 (cgs units) with 
10~ 2 solar metallicity. Estimates for the WD mass 
and radius were found by using theoretical cool- 
ing tracks on a T e fj - log(g) plane (Driebe et 
al. 1998). Figure 6 shows that the best fit T e ff 
and log(g) values give M w d — 0.47 ± 0.13M Q and 
Rwd — O.O165i? . Note that this is only a mass 
estimate since the WDs in CVs are heated by ac- 
cretion. The WD in PQ And might not follow this 
exactly. 

There was residual flux red-ward of 6000A be- 
tween the observed and the best fit WD model 
spectra. This red excess is from the secondary. 
To improve the model fit, a synthetic spectrum 
of an M dwarf was added to the best-fitting WD 
synthetic spectrum. Since the secondary in CVs 
fill their Roche lobe, their characteristics are de- 
termined by the properties of the binary system. 
Equation 2.101 in Warner (1995), 

i? 2 (0) = 0.094P o V° 8 (/i), (3) 

relates the orbital period and secondary radius. 
Applying this equation to our derived periods 
gives a range of R 2 = 0.17 - 0.20^©. The best fit 
effective temperature of the secondary from visual 
inspection was 2500 K, however, this should not 
be taken as definitive since the uncertainty is eas- 
ily ± 500 K. Even this range and the short orbital 
period still imply a late M dwarf as the secondary 
star. The best fit WD plus secondary model spec- 
tral energy distribution is shown in Figure 1. An 
attempt was made to extract the secondary spec- 
trum by subtracting the best-fitting WD synthetic 
spectrum from the observed spectrum. Unfortu- 
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nately, the resulting spectrum was too noisy and 
did not cover enough of the TiO and VO bands 
to make a reliable determination of the spectral 
type. The only reported infrared color for PQ 
And, (J - K) = 0.85±0.49 (Sproats, Howell, & 
Mason 1996), is consistent within the uncertainty 
with a late M spectral type. PQ And is too faint 
have been detected by the 2MASS survey (Hoard 
et al. 2002). 

The mass equation, Equation 2.100 in Warner 
(1995), 

M 2 (0) = 0.065P o 5 />), (4) 

is valid for P rb between 1.3 and 9 hours and pro- 
vides a secondary mass range of 0.13-0.16 M Q for 
periods between 1.7 and 2.1 hours. Given the 
derived WD mass the binary system mass ratio, 
q = M 2 /M u is 0.28-0.34. 

Using the derived WD radius, an absolute mag- 
nitude was computed for the best-fitting WD syn- 
thetic spectrum; M„ = 11.2. With an apparent 
magnitude, m.„ = 19 (Downes & Shara 1993) and 
E(B-V) = 0.06, the distance to PQ And was found 
to be 330 ± 50 pc. This distance places PQ And 
110 pc below the galactic plane. The only other 
distance estimate in the literature is from Sproats, 
Howell, & Mason (1996). They use the observed K 
band magnitude to determine the absolute magni- 
tudes. Their method depends on knowledge of the 
secondary's size and its fractional contribution to 
the K band flux. Without any period information 
Sproats, Howell, & Mason (1996) obtained a large 
range in distance for PQ And by assuming orbital 
periods between 80 minutes and 6 hours. Using 
their same method but now with our derived pe- 
riod range gives a distance between 319-817 pc. 
This estimate is still only an lower limit since it 
assumes that the secondary contributes 100% of 
the K band flux however it is consistent with the 
distance from our WD analysis. 

The best fit PHOENIX WD model puts PQ And 
in the region of the ZZ Ceti instability strip (Berg- 
eron et al. 1995). The low accretion rate implied 
by the strong Balmer absorption makes PQ And 
an ideal candidate for follow-up observations for 
non-radial pulsations. If found, an analysis of 
the oscillations will provide information on im- 
portant WD parameters including mass, compo- 
sition, magnetic field, rotation, temperature, and 
luminosity (Winget et al. 1991). GW Lib (van Zyl 



et al. 2004), SDSS J161033. 64-010223.3 (Woudt 
& Warner 2004), SDSSJ013132.39-090122.3 and 
SDSSJ220553.98+115553.7 (Warner & Woudt 
2004) are other WZ Sge type dwarf novae of simi- 
lar WD temperature that have already been shown 
to exhibit non-radial pulsations. 

5. Conclusions 

PQ And has characteristics similar to WZ Sge 
type variable stars. Its quiescent spectra show lit- 
tle indication of an accretion disk implying a low 
mass transfer rate. PQ And had three known su- 
peroutbursts of amplitude ~ 9 magnitudes in the 
50 years prior to 1988. Fits to spectra give P or b = 
1.7 hours while the period derived from photome- 
try is slightly longer at 2.1 hours. More observa- 
tions are required to confirm the true period but 
the available evidence strongly supports an orbital 
period below the orbital gap. The best fit WD 
model gives T eff = 12,000 K, log(g) = 7.7, and a 
distance of 330 pc. The WD parameters place it 
within the ZZ Cet instability strip. With a quies- 
cent magnitude of V = 19, non-radial pulsations 
could be detected in PQ And with a mid-sized 
telescope. 

SS is grateful to the NSF and NASA for partial 
support of this research. 
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Fig. 1. The combined spectrum (black line) and 
the best fit PHOENIX white dwarf model (grey line). 
The white dwarf model parameters are Teff = 
12,000 K, log(s) = 7.7, Z = 10~ 2 Z . The model 
also includes the flux from a M dwarf secondary 
with a T e// = 2500 K. 
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Fig. 2. — The evolution of Ha. The numbers on 
the right hand side of the figure give the time of 
mid exposure for each spectrum. NOTE TO ED- 
ITOR: This figure requires a minimum of 1/2 a 
page to show its details. 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 
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Fig. 3. — Radial velocity data obtained from a 
double Gaussian fitting to the line wings. Top 
panel gives the radial velocities as determined 
from the Ha emission feature. The bottom panel 
gives the radial velocities as determined from the 
H/3 emission feature. The best-fit solution yields 
a period of 1.7 hours. 
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Fig. 4. — The ratio of the blue (violet) peak to red 
Ha peaks. The solid line is the best sinusoidal fit 
assuming P or & = 1.7 h. 



0.40 E r 




1.10 F i .... i . , I , , , , 1 

197.65 197.70 197.75 197.80 

JD (-2451000.0) 

Fig. 5.— The 19 January 1999 differential light 
curve. The squares and diamonds are PQ And 
minus the comparison star and the comparison 
star minus the check star (plus 1 magnitude) val- 
ues plus their summed photometric errors, respec- 
tively. The line at the bottom left of the figure 
displays 1.7 hours in this plot. 
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Fig. 6. — The derived PQ And parameters on the 
theoretical white dwarf cooling tracks of Driebe et 
al. (1998). The White Dwarf mass is estimated to 
be « 0.5 M Q . 
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Table 1 
Emission line measurements 



Line 


EW 


Flux 


FWZI 




(A) 


(erg s _1 cm~ 2 ) 


(km s -1 ) 


m 


4 


2.6e-16 


2100 


H 7 


7 


4.6e-16 


2100 


H/3 


21 


1.0e-15 


2500 


He I 


5 


2.1e-16 




Ha 


75 


2.1e-15 


2700 



Table 2 

Radial Velocity measurements 







Ha 




H/3 


HJD a 


Phase 


RV 


Phase 


RV 






(km s _1 ) 




(km s _1 ) 


2898.890 


0.92 


149.1 


0.84 


116.6 


2898.897 


0.02 


-20.9 


0.94 


202.4 


2898.904 


0.12 


-23.4 


0.04 


72.8 


2898.911 


0.21 


-88.1 


0.13 


-86.4 


2898.918 


0.31 


-80.0 


0.23 


-88.2 


2898.925 






0.32 


-83.3 


2898.932 


0.50 


52.3 






2898.945 


0.68 


155.9 


0.60 


183.9 


2898.952 






0.70 


240.7 


2898.958 


0.88 


142.7 


0.80 


88.3 



a Hcliocentric JD of mid-integration minus 2,450,000. 
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Table 3 
PQ And system parameters 



Parameter Value 



Observed parameters 



RA (J2000) 


+02:29:29.54 


DEC (J2000) 


+40:02:39.40 


b 


-19.06° 


E(B-V) 


0.06 (Schlegel, Finkbeiner, & Davis 1998) 


Outbursts 


23 August 1938 (Richter 1990) 




7 March 1967 (Richter 1990) 




21 March 1988 (Hurst et al. 1988a) 


Recurrence Timcscalc 


25+4 years 


Maximum magnitude 


10 (visual) (Hurst ct al. 1988a) 


Minimum magnitude 


19.0 (V band) (Downes & Shara 1993) 


Ha emission line analysis 


Orbital period 


1.7+0.1 hours 


Systemic velocity 


46.0+3.5 km s" 1 (Equ. 1) 


Scmiamplitude 


135.7+13.7 km s" 1 (Equ. 1) 


Superior conjunction 


2,452,898.902 HJD (Equ. 2) 


A 


0.95+0.01 (Equ. 2) 


B 


0.25+0.01 (Equ. 2) 


PHOENIX WD analysis 


T eff 


12,000+1,000 K 


log(g) 


7.7+0.3 


Z 


10- 2 Z Q 


Distance 


330+50 pc 


Mi 


0.47+0.13 M 
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